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Quantum Mechanics SUper-SYmmetry (QM-SUSY) provides a general framework for studies 
using phenomenological potentials for nucleons (or clusters) interacting with a core. The SUSY 
potentials result from the transformation of the mean field potential in order to account for the Pauli 
blocking of the core orbitals. In this article, we discuss how other potentials (like external probes or 
residual interactions between the valence nucleons) are affected by the SUSY transformation. We 
illustrate how the SUSY transformations induce off-diagonal terms in coordinate space that play 
the essential role on the induced transition probabilities on two examples: the electric operators 
and Gaussian external fields. We show that excitation operators, doorway states, strength and sum 
rules are modified. 
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I. INTRODUCTION 

Almost all branches of many-body physics have devel- 
oped methods to simplify a many-body self-interacting 
system into a local "effective" mean potential affecting 
the pertinent degrees of freedom. It often provides an ad- 
equate starting point for more sophisticated approaches. 
For example, phenomenological potentials replacing the 
Schrodinger equation of N self-interacting particles by 
a one body potential whose orbits simulate the experi- 
mentally known structure have been widely used in nu- 
clear physics. Of particular interest are the halo systems 
which are often described in terms of valence nucleons 
interacting with a core. An elegant way to justify the 
effective core-nucleon phenomenological interaction is to 
invoke a super-symmetric (SUSY) transformation of the 
mean- field potential of a N-body system 0, S S 0, 0. 
Indeed, since the core is made of nucleons occupying the 
lowest orbitals of the mean-field potential, the halo nu- 
cleons cannot fill in these occupied states because of the 
Pauli exclusion principle. SUSY transformations includ- 
ing the forbidden states removal (States Removal Po- 
tential, SRP) and the restoration of phase shifts (Phase 
Equivalent Potential, PEP) provide an exact way to re- 
move the states occupied by the core without altering 
the remaining states properties. Hence, SUSY transfor- 
mation which can be fully analytical for some classes of 
potentials @ , provides an equivalent effective interaction 
between composite systems and thus can be safely used 
to describe nuclear structure and reaction of nuclei pre- 
senting a high degree of clusterization. 

SUSY transformations have been applied to breakup 
mechanisms involving halo nuclei |E0]- Indeed, the phe- 
nomenological treatment of halo nuclei in terms of nucle- 
ons in interaction with a core should take into account the 
fact that some intrinsic bound states of the nucleon-core 
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potential are Pauli blocked. For instance, the Is orbital is 
generally occupied by the core nucleons. In the case of a 
one neutron halo, like n Be or 19 C, SUSY-PEP potentials 
have been used to calculate -B(El) matrix elements 0], 
Coulomb breakup @], transfer reactions 0|. In the case 
of two neutrons halo, like 6 He, 11 Li, or 14 Be, it has been 
applied to remove the forbidden states and to analyze 
binding energies and radii of these nuclei [9. llfl llll 
Finally, it has also been included in coupled-channel cal- 
culations 0. El El- In an these calculations, SUSY 
transformation has been applied considering the follow- 
ing approximation: only the internal part of the Hamilto- 
nian (the core-halo potential) is SUSY transformed while 
the additional fields (external potentials, two-body cor- 
relations in the halo) remain unmodified. This approxi- 
mation will be called the internal SUSY approximation 
because it concerns only the core-halo potential. In the 
framework of this approximation, the SUSY transforma- 
tion is not equivalent to the exact treatment in which 
the Pauli blocked states are projected out. In this ar- 
ticle, we discuss a consistent SUSY framework which is 
always equivalent to the full projection-method. 

The accuracy of the internal SUSY transformation has 
been discussed in several papers. For instance, Ridikas 
et al. 0| have analyzed the radii of several halo nuclei 
as well as -B(El) matrix elements before and after the 
SUSY transformation. Thompson et al. ^3] an d Descou- 
vemont et al. 0] have performed a comparison of the full 
projector-method with the internal SUSY. As the consis- 
tent SUSY framework we discuss is totally equivalent to 
the full projection-method, the comparison between the 
internal approximation and the consistent SUSY treat- 
ment is thus an alternative method to estimate the accu- 
racy of the internal approximation. From the theoretical 
point of view, the consistent SUSY approach provides a 
unique and exact framework to compute excitation pro- 
cesses or to take into account residual interaction between 
valence nuclei. Such a consistent framework is essential 
to inter pre t the results of inverse problems in scattering 
theory [if}. 

This article is organized in the following way. In sec- 
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tion [HJ we develop a consistent formalism to map the 
original Hamiltonian into the SUSY partner Hamilto- 
nian. In the case of a static problem, this mapping is 
the usual one, but we will show in section UTA that in the 
case of an Hamiltonian modified by either an external 
field or a two-body interaction (for instance 2 neutrons in 
the halo), one should transform these fields into the new 
space. In section HVl we will illustrate the SUSY trans- 
formation showing both analytical results and numerical 
implementations for two potentials which are important 
in nuclear physics. We will then discuss the transforma- 
tion of external field: in section the response to an 
electric excitation of the general form ? x Ylm, and the 
response to a Gaussian potential in section IvTl 



II. SUSY TRANSFORMATION FOR ONE BODY 
HAMILTONIAN 



The application of smersymmetry to Schrodinger 
quantum mechanics [E 0- 0- Q- has shed new light 
on the problem of constructing phase-equivalent poten- 
tials. In this section, we review the SUSY transforma- 
tions which remove a state (SRP ) an d impose that the 
phase shifts are conserved (PEP) 0,0]. We will intro- 
duce mapping operators which change the Hamiltonian 
as well as the bound states. Finally we will present ef- 
fects on the external potential and residual interaction of 
composite systems described through effective Hamilto- 
nians. 



A. Initial Hamiltonian: ho 

Quantum mechanics SUSY has been extensively stud- 
ied for one dimensional systems. There are two ways to 
perform the multi-dimensional generalization depending 
on the choice of space coordinates. In three dimensions, 
one can choose the Cartesian coordinates, (x, y, z) [2( or 
the spherical coordinates, (r,0) We choose the lat- 
ter which is often used for excitation processes. Hence, 
the representation of the one particle Hilbert space TL is 
given by a sum over the sub-spaces 7i l associated to the 
angular momentum I: H — H° + H 1 + H 2 + . . . 

Let us introduce the initial Hamiltonian hg: 



ha = 

2m 



v , 



(1) 



where p is the momentum operator and where the po- 
tential operator v is assumed to be local. Since ho 
is rotationally invariant we can introduce angular mo- 
mentum as good quantum number and thus the wave 
functions associated to an energy E can be written as 
\<fa{E)) = i \(p l (E)) ® \Y hn ). In the sub-space H l , the 
radial static Schrodinger equation associated to the Zth 
partial wave is 



where p = ihV r is the radial momentum operator. The 
effective radial potential v l includes the centrifugal force 



h 2 1(1 + 1) 
2m f 2 



v . 



(3) 



In order to simplify the discussion, we do not include 
the spin-orbit potential. Nevertheless, the generalization 
of this framework to include spin-orbit potential is not 
difficult. 

In order to simplify the notations when there is no am- 
biguity we will drop the label I since the SUSY transfor- 
mations considered are defined in a subspace of angular 
momentum I (and m) i.e. they are block-diagonal in the 
complete space. Thus, they affect differently the poten- 
tial v associated with different I, 



B. Hamiltonian after k SUSY transformations: hu 

The elementary SUSY transformations remove a single 
state with or without restoring the phase shifts. In order 
to remove several states we will iterate the SUSY trans- 
formation. Therefore, let us assume that after k transfor- 
mations the radial static Schrodinger equation associated 
to the Ith partial wave can be written as 



P 

2m 



\ V \{E)) = E\v{{E)) . (4) 



It should be noticed that, since the SUSY transforma- 
tions are block-diagonal and different in each subspace of 
angular momentum Z, the different radial potentials do 
not correspond to the same potential i.e. the various v l k 
are different 



h 2 1(1 + 1) 



2m 



(5) 



We call E ( k l) (i = n, I) 
which is thus 



for different angular momenta I 
the energy of the ith bound state of h k 
(21 + l)-fold degenerate. 

The bound states correspond to the square integrable 
solutions of the differential equation (Eq. BJ. However, 
we will not restrict the solution of Eq.QJto bound states 
but rather consider all solutions \ifk(E)). Given a par- 
ticular solution \tpk(E)) of Eq. 0] whose inverse is square 
integrable, the general solution of Eq. 31 can be recast 
as 3] 



tpk(E,a;r) = ip k (E;r) 1 + a 



dr' 



, (6) 



where the parameter a can vary freely to construct all the 
possible solutions of the second-order differential equa- 
tion (Eq.QJ. In Eq.El we use the r- representation and the 
Dirac notations: £pk(E;r) — (r\(fk(E)). For future use 



^(E))=i—+v l A\^(E))=E\^(E)), (2) let us define the constant p = f J™ dr/& k (E;r)} 
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The Hamiltonians hk can always be factorized 

h = a k h a k + £ k , (7) 

where £k is the factorization energy and the first-order 

differential operators {aZ — (fit) ) are of the follow- 
ing form: 



'2m 



(hw k =F ip) 



(8) 



where w k = Wk(r) is the super-potential. Notice that 
in the literature, capital letters are usually used for the 
differential operators a h . Here, we dedicate capital let- 
ters for many-body operators while lower case letters are 
used for one-body operator. It is possible to show that 
the general solution \<p k (£ k, a)) of Eq. ^ with E = £ k 
is equivalently the solution of the first order differential 
equation 



a k \ip k (E 



■£k,a)) =0. 



(9) 



As a consequence, the super-potential is the local opera- 
tor defined by 

w k (E = £ k ,a;r) = — \nip k {E = £ k ,a;r). (10) 
ar 

For a given factorization energy £ k , there is a family 
of solutions which depends on the parameter a gener- 
ating the super-potential Wk(E — £ kl oi). Note that 
(pk(E = £k,a;r) must be nodeless in order af to be 
bound. Hence £k must be less than or equal to the ground 
state energy E k of hk and this requires also that a > —(3. 
The choice of the factorization energy £ k and the selec- 
tion of a member from the family of solutions w k must 
clearly be physically motivated. 

In this section we have defined the notations used in 
the following. In the next section we will present a 2 
step method which removes the lowest-energy state and 
preserves the phase shifts. 



C. State Removal Potential (SRP): v k +i 

The SRP transformation is defined so that it removes 
the lowest-energy state of a given sub-space H l . For the 
given angular momentum I, we choose £ k — E k , the en- 
ergy of the lowest energy state of the Hamiltonian h k . It 
follows that the inverse of the particular solution is not 
square integrable and it imposes a=0. With these defini- 
tions, we associate to hk a supersymmetric partner hk+i 
defined by 



K K ■+ 



p 

2m 



■ Vk+l, 



Vk 



(d r w k (£ k ,a = 0)) . 



(11) 
(12) 



The Hamiltonians h k and h k +x share the same spectrum 
except for the lowest-energy state of h k which has been 



suppressed in h k+ i. The states (\tp k +i(E))) of h k +i can 
be obtained from those (\ip k (E))) of h k according to 



\Vk+i(E)) 



1 



h k +i — £ k 



-.a k \<Pk(E)) = u k \>fik{E)) 



(13) 

Conversely, except for the ground state, the states of h k 
can be obtained from those of hk+i by 



\ME)) 



1 



\Jhk - £k 



%t\iPk + i(E))=n+\ i pk+i(E)). 



In these equations, we have introduced the pseudo uni- 
tary SRP-operators u^~ and u^" which are defined as (the 
products a k a k and a k a k being definite positive) 



ui = a. 



u fe = a k 



a k' a t V a fea fc 



k ■ 



(14) 



l k a k \ a k a k 



(15) 



These operators are pseudo-unitary since u k = (ii k )\ 

u fe u fe = 1 an d u fc u fc = Pi where- the projector p sup- 
presses the lowest-energy state \ip k ) of the Hamiltonian 
hk from the sub-space TC l and can be written as p = 

The relation between h k and h k+1 is 



hk 



fc+i 



P 

2m 



P 

2m 



Vk 



u+ h k u k 



(16) 

However, it is important to remark that ik+i ^ u^"wfcU fe 
and p 2 7^ u k p 2 u k . In fact the SUSY transformation of 
a local potential is not local. The simple diagonal form 
of the potential Eq.Elis recovered because, by construc- 
tion, the modifications of the kinetic part just cancel the 
off-diagonal terms in the transformed potential. Then, 
the kinetic and the potential parts of the Hamiltonian 
should be transformed together in order to get the rela- 
tion llKit with a simple potential (local in the r-space) 
and kinetic (diagonal in the ^-representation) terms. 



D. Phase Equivalent Potential (PEP): v 



■k+l 



It can be shown that SRP transformations change the 
phase shifts. To solve this problem, Baye Q has proposed 
to perform a second SUSY transformation and associate 
to h k+ i a new supersymmetric partner h k+1 so that 



A 2 
P 



hk+1 - a k a k +£_k— — 



v 



.fc+ii ' 



(17) 



with £_ k = E k , the ground state energy of h k> and 
a = —j3. In this case, the solution l P k+1 of h k+1 and 
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its inverse are not square integrable. The second SUSY 
transformation does not suppress nor add any state to 
the spectrum of h k+1 , but it restores the phase shifts so 
that the Hamiltonian h k+1 is equivalent to hk as far as 
the scattering properties are concerned. Note that the 
energy £_ k = E k is now below the ground state energy 

E%+i of hk+i. 

The corresponding super-potential w k (£_ k ) is deduced 

from the ground state wave function \(f1) of hk+i accord- 
ing to Eq. It can also be deduced from the ground 
state of hk according to 0] : 



dr' (/ fc (r')r 



m k (S. k ;r) - w k (£ k ;r), 



(18) 



where we have used the relation £_ k = £ k and introduced 
the modified-super-potential w_ k (£_ k ; r) as 

m k (£ k ;r) = ±ln £ dr> (rf(r')) 2 . (19) 
The corresponding potential v k+1 is 



Tr 



Vk+l 



(d r w k (£ k )) 



Vk 



hf_ 

rn 



(d r ± k (JLk) 



According to the above discussion, the spectra of hk+i 
and h k+1 are identical. All the states of h k , except its 
lowest-energy state, are mapped onto the states of h k+1 
with the same phase shifts. This mapping is simply: 



^+l( E ) 



1 



E° k -h k +i 



a k a k \tp k (E)) =u k \(f k (E)), 



E° k 



h k 



^+l( E ) 



(21) 
(22) 



with the pseudo-unitary PEP-operators 

1 1 



iiA- 



l k &k 



&k H 



a k a k 



(23) 



(24) 



In the following, as long as no confusion is possible, we 
will write the relations fulfilled by the general operator 
u k , which can be replaced either by the operator u k for 
the SRP transformation or ut for the PEP one. 



III. SUSY TRANSFORMATION FOR GENERAL 
HAMILTONIANS 

In nuclear physics, we are often interested in the de- 
scription of A interacting nucleons assuming that these 
nucleons can be separated into a frozen core contain- 
ing A c nucleons and a valence space containing A v nu- 
cleons. Hence, the wave function of this system is 
assumed to factorize into a core and a valence part, 
|$(A C + A V )) = |$c(A c )) ® |$v(A v )). The core state 
is described at the mean field level as a Slater deter- 
minant of A c single particle states \<j> c (ri)} occupying the 
A c lowest-energy eigenstates of the mean field potential 

ho: |$ C (A C )) = n«=i \ < f' c ( n )) w h ere ~ stands for the an- 
tisymmetrization sign. As a consequence of the Pauli 
principle, the valence nucleons cannot occupy the lowest 
orbitals of the core-valence potential which are already 
occupied by the core nucleons. The evolution of |$ V (A V )) 
is thus ruled by the Hamiltonian H v which contains a 
projection out of the occupied space P v = n«=i &n Ki 
(20) where c+ (c„) is the creation (annihilation) operator of 
a nucleon in the occupied orbital \<p c (n)). H v is assumed 
to contain the confining effect of the mean field ho- For 
cases with several nucleons in the valence space, the resid- 
ual interaction among valence nucleons, Vq should be 
taken into account when the problem of correlations is 
addressed. Finally, an external field, /o, should be intro- 
duced in order to compute excitation properties. Then 
the Hamiltonian reads 



PvffPv 



with 



A v A v A v 

H=j2h (i) + - £ v (t,j)+x;/o(i). 



(26) 



(27) 



i=i 



In the following, we propose to generalize the SUSY 
transformation so that it remains totally equivalent to 
the projector method for every kind of additional poten- 
tials. The first step of this method is to remove the k 
orbitals occupied by the core nucleons. We introduce the 
full operator t]^_ 1 which is the product of the different 
transformations u k _ 1 (c.f. Eq. ^Jand Ea. 12311 



The relation between h k and h k+1 is 



hk+i = U h k R k 



(25) 



The advantage of using the operators u k and u k is 
that all the relations we will deduce hereafter will be al- 
gebraically equivalent for SRP and PEP transformations. 



(28) 



applying to each single particle wave function (i) of the 
valence state. Since those different transformations af- 
fects only a given single-particle angular-momentum sub- 
space, the total operator U fc is block diagonal in spin 
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representation. Being the product of pseudo-unitary 
transformations, U fc is also pseudo-unitary since Ur = 



(U+j , U,TU+ = 1 and U+U- = P v . Using U+U~ = 
P v , we can thus write H v — P V HP V and explicitly 

= U+_ 1 U^_ 1 ffU+_ 1 U^_ 1 = U+^i^U^ (29) 

where we have introduced the transformed Hamiltonian 

A v A v A v 

= EM ! )+oE V k (iJ)+J^fk(i). (30) 



i=l 



It is clear from this relation that not only ho is trans- 
formed but also the two body interaction is changed into 
V k (i,j). Using Eq.Eland (i) u£ (i) = 1 (i) we get 

Mhj) = ^k-MiJ)^k-i (31) 

= K-i (*) u fc-i (j)Vo(h i) u fc-i 00 u fc-i (i) 

(32) 

and the external potential fo is mapped into 

= Ufe-i/o(*)Uti = ViW/oBfiliW. (33) 

Note that when /o is not a scalar operator, the mapping 
operator u^L^ on the right and on the left of Eq.EHlmay 
not correspond to the same angular momentum I. 

It should also be noticed that not only the Hamiltonian 
is changed but also the wave functions since the state of 
the valence nucleons is transformed into 



|*v,fc (*)> = U^ x |* v (*)) 



(34) 



The evolution of a state |$ v (*)) is driven by the time 
dependent Schrodinger equation 

\$y(t)) = H v |* v (i)) , 

which can be mapped into the new Hilbert space where 
the Pauli-blocked states have been removed within SUSY 
transformations: 

*n^|*y,fc(t)) = -Hy,* |*v,fc(*)) ■ 

It is important to remark that the projectors P v in- 
volved in the definition of the valence Hamiltonian (cf 
Eq. |22J have been removed in the mapped Hamiltonian 
H v ,k (cf En. l30j) . Hence, the time dependent Schrodinger 
equation in the SUSY space is simpler than the origi- 
nal Schrodinger equation which involves projection oper- 
ators. Nevertheless, the two Schrodinger equations writ- 
ten in the original space or in the SUSY transformed 
Hilbert space contain strictly the same physical ingredi- 
ents and are mathematically equivalent. 



In the literature, the transformation of both the exci- 
tation operators and the wave functions are usually ne- 
glected i.e. fo is often used instead of fk and the wave 
functions are not transformed back when evaluating ob- 
servables 0,0,I1I1E3E1E1E1E1 In the follow- 
ing, we will study a particularly important application 
which is the evaluation of the response of the nucleus to 
an external (one body) perturbation (time dependent or 
not) . The use of a SUSY transformed two body residual- 
interaction in the calculation of correlations and reactions 
will be the subject of forthcoming studies. 



IV. EXAMPLES OF SUSY 
TRANSFORMATIONS 

In this section, we illustrate the formalism developed 
above with two important physical examples: i) the har- 
monic oscillator potential which is mostly analytical, it 
allows a deeper insight into the formalism and provides 
numerical tests, and ii) the halo nuclei potential which 
are of important physical interest but can be treated only 
numerically since only asymptotic relations can be de- 
duced analytically. 



A. The harmonic oscillator potential 

The harmonic oscillator potential is a textbook ex- 
ample [ijj. We set the local potential to be Vo(r) — 
— Vq + ^( r /b) 2 with b 2 = h/muj. In the following, we 
introduce a reduced coordinate q = r/b. 

The eigenstates are labeled with the quantum numbers 
(n, I) and are associated with a set of energies E n \ = 
—Vq + (2n + I + 3/2)hto. For each I the lowest energy 
state is 

<p l (q) = ciq l+1 cxpj-iq 2 }, 

With Q = b 1 - 1 / 2 /^ 4 . 

We deduce the super-potentials for SRP and PEP 
transformations: 



w l (q) 



i fi+ 1 



b V q 
1 q 2l+2 erf 
b erf {q, 21 + 2) ' 



(35) 
(36) 



with the unnormalized error function defined by 
erf(z,l) = Jgt'e^dt. The differential operators a l k are 

(rloSV) = -( r \±(i±ip-hu l -±±)\r'), (37) 



V2 



(r& |r') S _ (r |_(_r±ip)|/) 



(38) 



where f = ^/mur and p = p/^/m. Notice that PEP 
transformation is not physical in this case because there 
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are no phase for the harmonic oscillator potential. Nev- 
ertheless, it remains interesting for the discussion. 

From the expressions of the super-potentials removing 
only one state, we deduce the transformed potentials: 




1 2 3 4 5 6 7 

Radius r [fm] 



FIG. 1: (color online) Radial part of the Harmonic Oscillator 
potential (H-O) for 1=0,1,2 (thick lines) compared with its 
SRP (dotted line) and PEP (dashed line) transformation of 
the 1=0 potential. The thin solid line stands for vu-o(i = 

l) + hw. 

These potentials are represented in Fig. In the 
graphical illustrations we will use nuclear physics scales 
by taking the following parameters: Vb=50 MeV and 
huj=10 MeV. The lowest energy state is at -35 MeV. The 
r.h.s. of Ea. 1391 demonstrates that the SRP transforma- 
tions removing only one state have mapped the original 
potential v l Q into a new potential which is simply v l + hcu 
where the effective angular momentum is I' = I + 1. This 
is illustrated in Fig. where we have represented the 
original potential with I = 0, I = 1 and I — 2 (thick 
lines) and the SRP potential obtained numerically (dot- 
ted line) . These numerical results have been obtained on 
a mesh containing 400 points, ranging from fm to 20 fm 
and with a vanishing boundary condition. The thin solid 
line is the analytical result given by Eq. EH1 The slight 
difference between the thin solid line and the dotted line 
gives an estimate of the error of the numerical algorithm 
which appears to be very small. 

Generalizing this result to the removal of several core 
states we remark that, in the new radial Hilbert space, 
up to a translation of n c hw, where n c is the number of 
removed core orbitals with the angular momentum I, the 
Pauli principle maps the original potential with angular 



momentum I to a new potential analogous to the radial 
potential with an effective angular momentum I' = I + 
n c . However, only the radial wave function is affected by 
the effective angular momentum, the angular part of the 
wave function is unchanged by the SUSY mapping. 

As we have already mentioned, this SRP transforma- 
tion changes the phases. The restoration of the phases 
is ensured by the PEP transformation. From the ana- 
lytic expression of (Eq. @HJ , we see that near r ~ , 

Uj (r) ~ Uq +2 ( r ) * an d asymptotically, v{ (r) ~ v l (r) . The 
potential v!i(r) is represented in Fig. (dashed line). The 
restoration of the phases imposes a non trivial transfor- 
mation of the potential: near zero, the potential v l is 
mapped to a new potential analogous to v l with a cen- 
trifugal force analogous to an effective angular momen- 
tum V = I + 2n c , and asymptotically, the potential stays 
unchanged as required by the phase conservation. This 
behavior is the consequence of the Pauli principle and 
phase restoration. 

The mapping operators, u^, can also be analytically 
derived, and we will discuss the properties of these op- 
erators from their asymptotic (all radii going to infinity) 
expressions: 

(rlufcV) £ -(r| = (ip + ?)[/) (41) 

y/P* /2m — -E/Q 

(HuoV) S S(r-r') (42) 

Hence, while the operator u^ is never trivial even at 
large distances the operator Uq* is simply the unity op- 
erator for large r. This is a consequence of the phase 
restoration. As a consequence, the PEP transformations 
do not modify observables which are only sensitive to the 
asymptotic part of the wave functions. These asymptotic 
properties are also valid for other potentials as illustrated 
for halo nuclei potential in the next paragraph. 

B. Halo nuclei potentials 

The study of the properties of weakly bound systems 
has found a renewed interest after the discovery of halo 
nuclei [20] . These systems have very large mean square 
radii and small separation energies. In fact, the separa- 
tion energy of the nucleons forming the halo is so small 
that their degrees of freedom can be separated from those 
of the nucleons constituents of the core. Up to now, this 
property has been observed in light nuclei close to the nu- 
cleon drip- lines like 6 He, n Be or 19 C. In reference HH, 
the proposed core-halo potential for 11 Be is the sum of a 
Wood-Saxon potential and a surface potential 

where f(r) is a Wood-Saxon potential and the parame- 
ters are: v = -44.1 MeV, a = -10.15 MeV, r\ = 1.27 
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fm and ciq = 0.75 fm. The bound states of this potential 
are 2 s-states at -25.0 MeV and -0.5 MeV and 1 p-state 
at -11 MeV. For simplicity we omit the spin-orbit cou- 
pling and consider a model case where the Is and lp 
orbitals are occupied by the core neutrons. Thus, these 
two orbitals are Pauli blocked and cannot be filled in by 
the neutron of the halo. In its ground state the latter 
occupies the 2s state. 



The creation/annihilation operators become 




2 3 4 5 

Radius r [fm] 

FIG. 2: (color online) Radial part of the core-halo potentials 
v l (solid lines) for I — 0, 1 and 2 and the SUSY transformed 
v[~° (dashed line) and v^~° (dotted line). 

We work on a constant step mesh containing 400 points 
and ranging from fm to 50 fm. We show in Fig. El 
the original potential v l for I = 0, I = 1 and I = 2 
(solid lines), the SRP v[=° (dashed line) and the PEP 
vj^° (dotted line). 

We can obtain analytical expressions near r — and 
for large r. Indeed, near zero, the lowest energy state 
behaves like r l+1 and asymptotically, it behaves like 
exp(— 7or) with 70 = \f— 2mJ^ /%. These asymptotics 
and therefore the following expressions are very general 
for all potentials which are regular at the origin and van- 
ish for large r. We find that the super-potentials behave 
like (r -> 0) 



~ and w l (r) ~ - 7o , (43) 

H^(r) ~ 2 ^ + 3 and jj^(r) ~ 0, (44) 



w (r) 
J 



, -j- 00 



. i 00 




(49) 
(50) 



Using these asymptotic expressions, one find the follow- 
ing properties of the mapping operators: 

(r\u \r) ~ (r| =M r ), ( 51 ) 



- v /p 2 /2m + 7^ 



<r|u±|r') SS 5(r-r'). 



(52) 
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FIG. 3: Matrix elements of Uq (r, r') (top, SRP) and Uq (r, r') 
(bottom, PEP) for several values of r' (r' = 0, r c , 2r c ) cal- 
culated for the angular momentum I = inside the core-halo 
potential. We remind that r c is the radius of the 10 Be core. 

We present in Fig. the matrix elements of (r|u^|r') 



and 



as a function of r for several values of 



and the potentials are for r — > 



0, r c and 2r c where r c is the radius of the 10 Be core. The 
peaks identifies the diagonal terms. We remark that the 
operator Ug~ has important off-diagonal terms (for small 
and large values of r') while the operator Uq converges 
towards a delta function when r' increases. Hence, the 
restoration of the phase shift imposes ~ 1 for large 
values of r', however this relation breaks down close to 
the core where the off-diagonal terms become important. 



ELECTRIC EXCITATION 



v\(r) A . (r) + ^^ = 4 +1 (r), 



to 

_2 ~ V \. r )i 



v[(r) Z v l (r) + 



and for 



v[(r) 
«i(r) 



v l Q (r). 



(45) 
(46) 

(47) 
(48) 



In this section, we shall consider dynamical properties 
of nuclei within the consistent SUSY transformation we 
have developed in the previous sections and discuss how 
it may be approximated. We will discuss the modifica- 
tions of the excitation operators, the doorway sensitivity, 
and finally, we will compute some transition elements and 
strength associated to monopole (E0), dipole (El) and 
quadrupole (E2) electromagnetic excitations and the as- 
sociated sum rules. 
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We assume that, prior to any SUSY transformation, 
the excitation operator takes the standard multipolar 
form 



/o(A,L,M) = / rad (A)r L , M , 



(53) 



with the radial excitation operator io ad (A) = ? X - We 
drop the coupling constant because we are only inter- 
ested in the transformation of the radial excitation op- 
erator and the relative difference between the consistent 
SUSY transformation and its approximations. The EO 
transition is induced by /o(2,0, 0) and the electromag- 
netic transitions EA are induced by /o(A, A, M) (A > 1). 
The SUSY transformation of the excitation operator is 

A(A, L, M) = u+_ x / (A, L, M) <£_ x . 

Introducing explicitly the angular momentum quan- 
tum numbers, the radial excitation operator is thus given 
by 



if(A) = u<'+ 



!/5 ad (A)ui-_ 



(54) 



where u^._ x is the mapping operators u fc _ 1 in a subspace 
associated with the angular momentum I (and to). The 
external operator (A) allows transitions between dif- 
ferent angular momentum space according to the selec- 
tion rules deduced from the relation 

(l'm'\f k (X,L,M)\lm) = /£' : (A)(/' 'm'\Y L>M \lm) . (55) 

It should be noticed that, in Eq.Q2|the mapping operator 
on the right and on the left side of /o ad (A) may not 
correspond to the same angular momentum I, Moreover, 
while the original radial excitation operator /($ ad (A) is 
diagonal (in the coordinate space), (A) is no longer 
diagonal because the transformation operators uj^ are 
non local. 



a SRP or a PEP transformation to remove this occupied 
state from the halo Hilbert space and restore the phase 
shift. Of course, to be complete we have also to remove 
the occupied p-state but since the SUSY transformation 
is block-diagonal for the angular momentum quantum 
numbers this does not modify the dynamics in the s- 
subspace. 
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FIG. 4: (color online) The normalized difference of diago- 
nal matrix elements (-R fc las (A; r) for SRP transformation and 
_R fc las (A; r) for PEP one) are represented for several values of 
A (1,2,4,8) for the electric excitation -EA inside the s-space. 

In order to evaluate the difference between the con- 
sistent SUSY transformation and its approximations, we 
define two quantities 



i?fs(A;r) = 
Rf(X;r,r') = 



yjA/°°(A)|r) 
(r|/HA)k) 

(r\f2°(X)\r) ' 



(56) 
(57) 



A. Consistent excitation operator and its 
approximations 

In the following, we shall perform the calculation of 
the excitation operator and some of the observables it 
induces. In the literature, the SUSY transformation is in 
general not applied to the excitation operator. Hence, in- 
stead of calculating the matrix elements induced by the 
consistent excitation operator f l k l (X), the authors have 
evaluated the matrix elements of /g ad (A) with the SUSY 
transformed wave functions. We will refer to this approx- 
imation as the internal approximation. We introduce a 
second approximation called the diagonal approximation 
which consists simply of neglecting the off-diagonal terms 
in coordinate space of the consistent excitation operator. 

As a first example we will study the -EO excitation of 
the halo neutron in the s-subspace. In this subspace the 
core blocks one orbital (Is) so that we have to perform 



where 



A/£*(A) = /£'(A)-/r(A) 



ft'i 



is the difference between the excitation operator consis- 
tently transformed f l k 1 (A) and the original excitation op- 
erator /Q ad (A). The ratio i? dlag (A;r) evaluates the dif- 
ference between the diagonal part of the consistent ex- 
citation operator and the original operator, normalized 
to the value of the diagonal part of the consistent oper- 
ator. It gives an evaluation of the approximation for the 
diagonal part of the excitation operator. Fig. 31 shows 
the ratio i? diag (A;r) and ^ ias (A;r) for A=l, 2, 4, 8 
in the case of the Is SRP and PEP transformation re- 
spectively. We remark that for large radii r, the diago- 
nal part of the excitation operator is close to the orig- 



inal one (i? fe lag (A; r) ~ 0). This is a consequence of 
the asymptotic properties of the mapping operator as 
it has been discussed in the section IIV Bl For small 
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radii r, the diagonal part of the excitation operator is 
strongly modified, the ratio J2 fc lag (A; r) ~ 1 revealing that 

( r l/fe ad | r ) ^ ( r l/o ad | r )- This difference persists through 
a large range of radial coordinates. This range increases 
with A, it is wider for SRP transformation compared with 
PEP transformations. 



r'=2r 



defined as 



FIG. 
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PEP -_ 





[fm] 
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5: The normalized off-diagonal matrix elements 
(7?fe tt (A = 4;r, r') and i?fc ff (A = 4;r,r')) are represented as 
a function of r for several values of r' (0, r c , 2r c ) for the 
electric excitation EX inside the s-space. 



\5MW^i')) = fkWk) 



In the following, we have chosen for \ip l k ) the ground state 



of hi, 



8(p.(2,0^0;r) 
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FIG. 6: (color online) The doorway state 5ipi(X = 2, — > 0; r) 
and 5^ (A — 2,0 — * 0; r) of electric excitations for SRP and 
PEP transformations respectively. The solid line stands for 
the consistent excitation operator, dotted line for the internal 
approximation and dashed line for the diagonal approxima- 
tion. 



On the other hand, the ratio i?£ ff (A; r, r') evaluates 
the amplitude of the off-diagonal terms in coordinate 
space normalized to the diagonal term of the consistent 
excitation operator. In Fig. we represent the ratio 
Rf(X = 4;r,r') and Rf(X = 4;r,r') for the SRP and 
PEP transformation respectively, and for three values of 
r': 0, r c and 2r c . For SRP transformation, off diagonal 
terms are important for small r' and decrease in relative 
magnitude while r' increases. Off-diagonal terms are non 
zero in a wide range and we will show in the next para- 
graphs that they can have a more important effect on 
observables than the diagonal terms. For the PEP trans- 
formation the off diagonal terms are smaller and become 
negligible for intermediate and large r' (> r c ) as required 
by the restoration of the asymptotic behavior. 

In all the cases presented here, the consistent excita- 
tion operator is different from the original one in the 
space region inside the core potential. Hence, from this 
observation, we can expect that there will be important 
effects induced by the consistent calculation if and only if 
the calculated observable is sensitive to the space region 
inside the core potential. 



B. Transformation of the doorway state 

We want now to evaluate both contributions of the di- 
agonal and off-diagonal excitation operators on the ma- 
trix elements. For that, we introduce the doorway state 



We represent in Fig. the doorway state 6<pi(X — 
2,0 -► 0;r) and dip^X = 2,0 -> 0;r) for the SRP and 
PEP transformations respectively. The solid line stands 
for the consistent excitation operator, dotted line for the 
internal approximation and dashed line for the diagonal 
approximation. We remark that the internal approxima- 
tion and the diagonal approximation are indistinguish- 
able. This shows that off-diagonal terms are the most 
important sources of modification of the excitation oper- 
ator. 

Moreover, the consistent doorway state changes sign 
while the two approximations remains positive. This af- 
fects the node structure of the wave function and may 
induce strong modifications for forbidden transition as 
we will see in the next paragraph. 

C. Single particle reduced transition probability 

The single-particle reduced transition probabilities are 
defined as [22] 

B k (E0,i^f) = ( V {\f l / h (2M) (58) 
with If = h and as 

B k (EX,i^f) = (<p{\f l k fh (\M) (59) 

with \lj, — A| < If < h + X and for A > 1. In order 
to simplify the notations, the states i and j are labeled 
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according to the original space prior to any transforma- 
tion. In the halo case, developed above, all the final 
states are in the continuum so it will not be possible to 
use directly these definitions. In the next section we will 
introduce the strength function, a more general way to 
look at transition probabilities which is suitable for the 
case of excitation towards continuum and which can thus 
be used in the halo case. To get results for the transi- 
tion probabilities between discreet states, in the present 
section, we will restrict the discussion to the harmonic 
oscillator model (see section llV All . To simplify the dis- 
cussion, we will consider that the nucleons of the core 
only occupy the Is orbital and we will study the exci- 
tation of an additional neutron in the 2s orbital. We 
have computed numerically the reduced matrix elements 
Bi(EO) and ^(El) for the PEP transformation. The re- 
sults are presented respectively in the tables U and [H] In 
the harmonic oscillator, due to selection rules, from the 
2s the monopole operator f 2 can induce transition only 
towards the 3s. In table HI the first line shows the result 
of the matrix elements (B) induced by the consistent ex- 
citation operator. As expected, the forbidden transition 
are zero within the numerical uncertainty indicated in 
parenthesis. 

The matrix elements, B, induced by the internal ex- 
citation operator are showed in the second line of table 
HI For allowed transitions, the internal approximation 
modifies by about 20% the exact matrix element, but 
the main effect of this approximation is that it induces 
forbidden transitions from 2s to 4s-7s states. 

On the other hand, in the case of El electromagnetic 
transition, the selection rules of dipole transitions in the 
harmonic oscillator allow transition from 2s states to lp 
and 2p states. The same phenomenon is observed in ta- 
ble U and table |nj the internal approximation produces 
spurious excitation of forbidden states. Moreover, for the 
allowed transition the error goes up to more then a factor 
2. 



f 


3s 


4s 


5s 


6s 


7s 


B^EO) 


l.lxlO 4 


o(10" 4 ) 


o(10- 8 


) o(10" 7 ) 


o(10" 7 ) 




9.2xl0 3 


1.5X10 1 


1.3 


7.1xl0" 2 


1.4xl0" 4 



TABLE I: B^EO, 2s-^f) for the PEP transformed harmonic 
oscillator potential. B^EO, 2s— + f) is the matrix element 
induced by the internal approximation of the excitation op- 
erator. 



D. Strength 

The very important discrepancy between the internal 
and the complete SUSY observed in the case of the har- 
monic oscillator might be a peculiarity due to the sym- 
metry of this model. Let us thus come back to the phys- 
ical case of the 11 Be halo nuclei for which transitions 
between orbitals belonging to the same 1-space are all al- 
lowed. Now we should remove the Is (-25 MeV) and lp 
(-12 MeV) orbitals occupied by the core neutrons so that 
only one bound state (2s) is available for the halo neu- 
tron. The excitations can only promote the halo neutron 
to the continuum. Hereafter, the eigenstates and the con- 
tinuum states will be obtained from the diagonalization 
of the Hamiltonian inside a box going up to 50 fm with 
400 points. 

In order to discuss transition towards the continuum, 
we introduce the strength: 



Sfc(EA,i, w) 



f)\ 2 6(LJ-E f +Ei), 



(60) 



where i is the initial state, here the 2s orbital, and the 
final states / are the continuum states of the box. The 
single particle energies are Ei and Ef respectively. Since 
we perform the calculation in a box the continuum is dis- 
creetized. In order to obtain a smooth strength function 
one often smoothes the obtained results with a Gaussian 
or a Lorentzian function. In this paper we will do both. 
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5.0xl0 2 1.2X10 3 
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ii(Ei) 


1.2xl0 3 8.0xl0 2 
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1.2 


7.8xl0" 2 



TABLE II: B^El, 2s^ f) for the PEP transformed harmonic 
oscillator potential. B^El, 2s^ f) is the matrix element in- 
duced by the internal approximation of the excitation opera- 
tor. 



FIG. 7: (color online) For monopole excitation: Left (a) 
the bars are the reduced transition probabilities Bj^EO) for 
the continuum states discreetized in the considered box, the 
solid line is the strength function resulting from a Lorentzian 
smoothing (F=500 keV) and the dashed line is the result of 
the internal approximation. In part (b), we show the ratio 
^(EOVl^EO). 

Strengths for PEP transformations for E0, El and E2 
transitions are presented in the parts (a) of Fig. I7I8I9I 
with a Lorentzian smoothing (r = 500 keV) and Fig. 
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FIG. 8: (color online) Same as Fig.|3for dipole transitions. 




Energy Co [MeV] Energy go [MeV] 

FIG. 9: (color online) Same as Fig. for quadrupole transi- 
tions. 

with a Gaussian smoothing. The part (b) of each figure 
gives directly the ratio B_ x (EA ) /B_ x (EA ) computed for 
individual states. 

For the monopole mode, using a Lorentzian smooth- 
ing the consistent strength and the internal strength ap- 
pears to be very similar (see Fig. 0-a) even if the ratio 
B_ x (EO) /B_-y (EO) computed for individual states (see part 
(b) of the figure) is very different from 1 for large values 
of the final excitation energy Ef. 

The dipole excitations connect the states of two differ- 
ent 1-subspace. On the smoothed S-^El) strength (see 
Fig. 0-a) we only observe a small over estimation of the 
strength for large values of Ef but again the effect seems 
much larger on the ratio B 1 (E1)/ ; B 1 (E1). In Fig. we 
represent the strength S_ ± (E2) . In the 1=2 subspace, there 
are no core states and consequently, the SUSY trans- 
formation is in fact the unity. The smoothed strength 
appears to be only slightly under-estimated by the inter- 
nal approximation again in contradiction with the ratio 




Energy co [MeV] Energy co [MeV] 

FIG. 10: (color online) Same as Fig.|Uwith G aussian smooth- 
ing. 

£? 1 (E2)/£? 1 (E2) which exhibits a strong discrepancy. 

In order to solve the contradiction we have first stud- 
ied the role of the smoothing. We have found that the 
situation is different with a Gaussian smoothing as illus- 
trated in Fig. This difference is due to the difference 
in the tails of the two smoothing functions: the long tails 
of the Lorentzian function associated with the low en- 
ergy states which have a large B_ x (EO) dominate even at 
large energy when a Lorentzian form factor is used. In- 
deed, since the difference between the two calculations 
is small for these dominating states the final Lorentzian- 
smoothed strengths for the two calculations are rather 
close in contradiction with the direct ratio of individ- 
ual excitation probabilities or the results of the Gaussian 
smoothing. 

To avoid the ambiguity of the smoothing method we 
have studied the continuum limit by a direct scaling of 
the numerical box size (we have also tested the role of the 
mesh size) . We have observed that the ratio B_ x /B^ com- 
puted for each individual state does not change shape go- 
ing to the continuum limit while the smoothed strengths 
vary and exhibit a strong dependence into the smooth- 
ing functional and parameters. Therefore, the ratio 
B_ 1 (E0)/B_ 1 (E0) provides in fact the correct continuum 
limit and the large observed discrepancy at high energy 
between the internal and the complete SUSY transforma- 
tion is the physical one. This is even better illustrated 
by considering integrated effects like effects on the sum 
rules. 



E. Sum rules 

By integrating the strength over the energy, we can 
define different sum rules m t 

m t (EX) = [du ui^iEX, i = 0, u), (61) 
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where t is the weight of the energy. In the frozen core 
approximation, the sum rules mo and mi can be obtained 
from the halo wave function according to 

m = ($v|/oPv/o|$v) - (<$ v |/o|$v)) 2 , (62) 
mi = i($ v |[/ Pv,[^v,/oPv]]|$v>. (63) 

where fo is the excitation operator defined by Eq.E3 We 
define the internal approximation for the sum rule as fht 
for which the projector P v have been removed. The Pauli 
principle is no longer respected. In order to estimate the 
error induced in the calculation of fh t compared to m*, 
we have estimated the ratios (m t — m t )/m t and we have 
found the results presented in Table ITTT1 The relative er- 
ror induced by the internal approximation increases with 
the weight. This result is compatible with the results pre- 
sented in Fig. and Fig. when the weight increases, 
the contribution of large energy increases and the dis- 
crepancies between the consistent SUSY and its internal 
approximation increase also. 





(mo 


-pepm )/m (m^ 


-Mi) Mi 


(m 2 -m 2 )/m 2 


EO 




1.4% 


3.9% 


17.4% 


El 




-6.8% 


-33.3% 


-93.1% 



TABLE III: This table presents the relative error induced 
in the calculation of the energy weighted sum rules by the 
internal approximation of PEP transformations. 



VI. RESPONSE TO A GAUSSIAN EXCITATION 



as 

i x ffo ( (r- r ) 2 \ 

where the norm of the interaction is <7o=450 MeV.fm 3 
and its range is /i=2 fm. The SUSY transformation of 
this potential is 

Similarly to the previous section, we define two quanti- 
ties that measure the modification of the SUSY trans- 
formation on the diagonal and off-diagonal terms of the 
Gaussian potential: 




In the previous section, we have shown that the PEP 
transformation modifies essentially the external excita- 
tion operator in the space region located inside the core 
potential. The electric operators f x YLM studied in the 
previous section which can be seen as a multiple expan- 
sion of a Coulomb field far from the nucleus or as the 
low momentum transfer limit of a plane wave scattering 
are strong far from the nucleus. Hence, the effect of the 
PEP transformation on the excitation process has been 
found to not be too large. However, this is not always 
the case and in particular nuclear scattering and/or large 
momentum transfer reaction correspond to much shorter 
distances. In order to study the effect of the PEP trans- 
formation on such kind of scattering, in this section, we 
study the response to a Gaussian excitation which can 
strongly overlap with the core potential. Gaussian po- 
tential can be induced by an external nuclear potential 
as well as a residual two-body interaction between parti- 
cle in the halo. In a similar spirit of the previous section, 
we will not study a specific process but the response to a 
one body Gaussian potential centered around Tq defined 



_inLZ-j — , — , — , — I — i — i — i — , — i _n 5 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
5 10 5 10 

Radius r [fm] Radius r [fm] 

FIG. 11: (color online) Left (a) the initial Gaussian 
potential^ = 0) (solid line), SRP (dotted line) and PEP 
(dashed line) potentials. In part(b), we show the ratio 
Qf ag (r) for SRP (dashed) and PEP (dotted) potentials. 

In Fig. El we fix r o = and we represent (a) the di- 
agonal part of gi (dotted line) and g 1 (dashed line) com- 
pared to the original potential go (solid line), and (b) the 
ratio Q^ las (r) for the SRP transformation (dashed line) 
and PEP transformation (solid line). In the very cen- 
tral region, the PEP transformations modify the poten- 
tial by a factor 30%. At large distance r, because of the 
Gaussian shape centered on zero of go( r ')> an important 
relative weight is given to small radii in the summation 
( r \9k\r) — J dr'(r|u - |r')(7o(?"')( r/ | u l r ')- The result of this 
effect is that the range of (r\gk\r) is slightly increased 
compared to go(r') and because of the exponential be- 
havior of the excitation operator this is enough to make 
the ratio Qf &s (r) ~ 1. 
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FIG. 12: The ratios Q° e (r,r') in the I = channel as a 
function of r for several values of r' (0, r c and 2r c ) in the case 
of a Gaussian potential. 



In Fig. Q21 we show the ratio Q° s {r,r') for the SRP 
transformation (upper panel) and the PEP transforma- 
tion (lower panel) and for three different values of r'; 
0, r c and 2r c . For values of r' inside the potential, off- 
diagonal terms are small compared to the diagonal term 
but they are spread over a large range of coordinates, and 
their integrated effect can counter balance their small val- 
ues. Outside the potential, the off-diagonal terms become 
very important and even larger than the diagonal term 
for both SRP and PEP transformations. 

Both diagonal and off-diagonal terms have an effect on 
the particle wave function that can be estimated with the 
doorway state 

\5<p k {l^l'))=gl l \<p l k ). 

In the following, we have chosen for \ip l k ) the ground state 
of h[. 

In Fig. El we represent 6(fi(0 — > 0;r) for SRP and 
PEP transformations. The solid line stands for the con- 
sistent transformation of the Gaussian interaction, the 
dotted line for its internal approximation and the dashed 
line for the diagonal approximation. The internal ap- 
proximation and the diagonal approximation give about 
the same results but are both very different from the 
consistent calculation. These two figures illustrate the 
importance of the off-diagonal terms which induce very 
different doorway states. 

Summarizing our results, we can assert that the Gaus- 
sian interaction is considerably modified by the SUSY 
transformations and the internal approximation is cer- 
tainly a bad approximation as it is illustrated in Fig.lT3l 
Hence, calculations of structure properties and reaction 
mechanism which involve SUSY transformations should 
never neglect the transformation of the excitation op- 
erator (or residual interaction) for radii inside the core 
potential. 
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FIG. 13: (color online) The doorway state in the I = channel 
8(fii(0 — > 0; r) produced by a Gaussian excitation for SRP and 
PEP transformations. The solid line stands for the consistent 
excitation operator, dotted line for the internal approximation 
and dashed line for the diagonal approximation. The thin line 
stands for the initial wave function ip° (r) . 



VII. CONCLUSION 

In this article, we discussed a consistent framework 
to perform quantum mechanics SUSY transformation to 
take into account the internal degrees of freedom in a 
core approximation. This method is totally equivalent to 
the full projector-method and is formally very interesting 
since it provides justifications for effective core-nucleon 
interactions and nucleon-nucleon residual interaction. In 
this study, we have considered several kinds of external 
fields and performed a consistent SUSY transformation. 
The consistent SUSY transformation provides equivalent 
effective interactions between composite particle systems 
and thus can be safely used to describe nuclear structure 
and reaction of nuclei. The consistent transformation of 
additional fields as well as the transformation of the wave 
functions (or the observables) is usually neglected in the 
literature (internal approximation) and we have shown 
that it is not always justified. 

Our conclusions are the following: in the case of elec- 
tromagnetic induced transitions, consistent SUSY trans- 
formation conserves all selection rules while the inter- 
nal approximation violates it. Performing different com- 
parisons we have shown that the discrepancies might be 
large, affecting the node structure of the doorway states 
and changing the transition probabilities by sizeable fac- 
tors. Even the sum rules can be affected by a large per- 
centage, e.g. 33% for the energy weighted sum rule of 
the dipole excitation. Hence, the use of the internal ap- 
proximation for external excitation operator might be 
dangerous and the results obtained should be carefully 
discussed. But the main discrepancies between the con- 
sistent calculation and its internal approximation appear 
for external fields which have a strong overlap with the 
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core potential. For instance, it is the case of the Gaussian 
interaction centered at small distance (r < r c ). 

In all the cases, we have shown that due to the SUSY 
transformation, the off-diagonal terms of the external 
fields are often more important than the diagonal ones. 
This forbids an approximation which would be to take 
into account only the SUSY modification of the diagonal 
term. Hence, the SUSY transformation has to be fully 
implemented in order to preserve the symmetry of the 
original Hamiltonian. 

All the discussion related to the excitation operator is 
valid for the observables. Since the wave functions are 
transformed either they should be transformed back be- 
fore evaluating average values or the observables should 
be also transformed before being applied on a trans- 
formed state. 

In conclusion, in this article, we have stressed the im- 
portance to keep the consistency of the QM-SUSY frame- 
work if there is an overlap between the core potential and 
the additional interactions (excitation operator or ob- 
servables). For instance, in a recent article Hesse et al. |j| 



have performed the internal SUSY approximation and 
they have shown that in order to reproduce the known 
binding energies and radii of 6 He, 11 Li and 11 Be halo 
nuclei, a readjustment of the core-neutron interaction is 
required. This effect might be induced by the internal 
SUSY approximation which treats improperly the r 2 ob- 
servable for the halo neutrons and the residual interaction 
between them. The consistent SUSY framework would 
be a way to extract informations concerning the neutron- 
neutron interaction in the halo because there is a unique 
mapping between the original known interaction and the 
effective one which includes consistent removal of core 
orbits. 
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